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Abstract 

It is shown that the weak phase 7 = arg(—V u dV* b V c bV* d ) can be determined using only untagged 
decays B°/B° — ► DK$. In order to reduce the uncertainty in 7, we suggest combining information 
from B^ — ► DK^ and from untagged B° decays, where the D meson is observed in common decay 
modes. Theoretical assumptions, which may further reduce the statistical error, are also discussed. 
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I. INTRODUCTION 

CP violation measured in B — > J/ipKs |l|] is interpreted in terms of the phase (3 
axg(—VtbV t * d V c dV* b ) in a way which is practically free of theoretical uncertainties, providing 
an important test of the Kobayashi-M askawa mechanism Q . On the other hand, the current 
interpretation of CP asymmetry measurements in B — ► 7r + 7r~ 0] in terms of the phase 
a = &rg(—VtdV t * b V u bV* d ) involves an uncertainty in the ratio of penguin-to-tree amplitudes Q|. 
A theoretically clean method 0,0] for measuring the phase 7 = &rg(—V u dV* b V c bV* d ) involves 
interference between tree amplitudes b — > cus and b — > wcs, governing 5 — > D°X S and S — > 
D X s , where the _D° and _D° decay to a common hadronic state, and X s = K, K*, Kw, . . . 
is a strangeness one state. 

Originally, this idea for measuring 7 was proposed for charged B decays 5] and for 
time-dependent neutral B decays 6] of the type B — > D CP K. Here positive (negative) CP 
eigenstates, such as K + K~ (Kstt ), identify equal admixtures of D° and D° states with equal 
(opposite) signs. A variety of other D decay final states can be used as well, leading to several 
variants of the original method 0,0. Every hadronic state accessible at tree level to D° decay 
is also accessible at tree level to D° decay with varying levels of Cabibbo suppression. Flavor 
states, K~tt + and K*~tt + , which are Cabibbo-favored in D° decays, are doubly Cabibbo- 
suppressed in D° decays. Flavorless states, such as K* + K~ and K + K*~, are produced in 
singly Cabibbo-suppressed decays of both D° and D°. Recently, it was shown that a model- 
independent extraction of 7 is also possible by considering B^ — > DK^ with subsequent 
multibody D decay, such as D — ► K s tt + tt~ 0, liol ]. Other variants make use of the decays 
B± -> TfK^.D* -> Dn°, the self-tagged decay mode B° -> DK*°,K*° -»■ K + 7T~ Jli, 



multibody B decays of the type B — > DKir 12| and combinations of these processes 



First results for B^ — > Dq?K ± were presented recently by the Belle ^| and BaBar[n 
collaborations. These studies were based on several tens of events in each experiment and 
demonstrate the potential of a larger data sample in providing useful constraints on 7 jl^ . 
The Belle collaboration has also presented a preliminary analysis of about one hundred 
events of the type B ± -> Z)^, D -> K s n + 7T [17| . from which constraints on 7 were 
obtained. The main difficulty of each of these methods is that each decay mode by itself 
has a very low rate. Reaching high sensitivity in near future measurements of 7 requires 
combining several relevant B and D decay modes. 

Measurements relevant to studying 7 have so far focused on charged B decays, — ► 
DK ± and B ± — > DK*^ 1 At first glance, neutral -B decays seem to be less promising 
for two reasons: they have much smaller rates and require B° flavor tagging. Let us discuss 
these two points one at a time: 



The processes B° — > DK° and B° — > DK*° are expected to be color-suppressed 11], 
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implying that their rates are about an order of magnitude below the rates of corre- 
sponding charged B decays. However, the crucial factor determining the sensitivity 
of a measurement of 7 is not the decay rate itself. Rather, the sensitivity is governed 
by the magnitude of the smaller of the two interfering amplitudes in a given process. 
Since the smaller amplitudes in B + and B° decays are both color-suppressed, they are 
expected to be of comparable magnitudes. Therefore, this by itself is not a limiting 
factor for neutral B decays. 

• Only time- integrated rates have so far been measured in B° — > DK°, combining rates 
for B° -> D°K° and B° -> D°K° Q- It was shown in 0,0 El | that a determination 
of 7 is possible from time-dependent measurements, for which one must tag the flavor 
of the initial B°. The effective flavor tagging efficiency at B factories is about 30% 
(and much smaller at hadron machines), resulting in a doubling of the statistical error 
relative to the perfect-tag case. As we show below, 7 can be determined using untagged 
data alone. This makes use of events that cannot be tagged or even events that are 
mis-tagged, regaining a significant part of the sensitivity lost due to the low effective 
tagging efficiency 

In the present paper we investigate what can be learned from untagged decays B°/B° — > 
DKs, where the D meson is observed in several decay modes. A potential lower bound 
on I cos 7 1 from unt agg ed decays, in which D mesons are observed in CP-eigenstates, was 
noted by Fleischer 22]. We will go beyond this bound by showing that 7 can actually be 
completely determined in the range < 7 < n, using only untagged decays. In practice, 
it is useful to combine information from untagged neutral B decays with information from 
charged B decays, since the observables related to D decays are common to both cases. 
This provides an overconstrained information, permitting a more accurate determination of 
7 than when using B + decays alone. 

The plan of the paper is as follows. In Section II we introduce parallel notations for 
charged and neutral B — > DK decays, discussing briefly relative magnitudes of decay am- 
plitudes in these processes. Section III studies two-body and quasi two-body D decays, 
distinguishing between several classes of decay modes. We show that 7 can be determined 
from untagged neutral B decays alone and derive an explicit expression for tan 2 7 in terms of 
measurable rates. Multibody D decays in B° — > DKs are studied in Section IV. In Section 
V we discuss a way of reducing the number of hadronic parameters by assuming isospin 
symmetry and by neglecting an annihilation contribution in B — > DK. Finally, Section VI 
concludes. We also add an appendix studying time-dependence in B°(t) — > (K s i{ + i{~)dKs- 
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II. AMPLITUDES IN B —> DK DECAYS 



We define decay amplitudes of B — > DK for charged B mesons, 



A(B + -> D°K + ) = A c , 

A(B + -> D°K + ) = A c r c e i{5c+ ^ , (1) 



and for neutral B mesons, 



A(B° - D°K S ) = A n , 

A(B° - D°K S ) = A n r n e 1 ^ . (2) 

By convention, Aj > 0, > 0, and < 5i < 2tt (i = c, n). Amplitudes for the CP conjugated 
decays have the same expressions, but the phase 7 occurs with an opposite sign. 

Let us discuss briefly the relevant ratios of amplitudes. The amplitude A c is a combination 
of color-allowed and color-suppressed contributions, while the amplitude A n is purely color- 
suppressed. The ratio A n /A c may be estimated in two ways leading to comparable values. 
We mention in each case the required approximation: 

• Measurements of B + D°K+ Q, H Q and B° -> D°K° imply 



In addition to a term A 2 nl the untagged rate in the numerator includes also a smaller 
term A 2 n r\ which we neglect. 

Using flavor SU(3) [24], one may relate A n /A c to a corresponding ratio measured in 

B ->• Dir Ei, 

' V -0.24 ±0.01. (4) 



This relation is affected by SU(3) breaking corrections, and by a small exchange am- 
plitude in B° — >• D°tt° [26, 27]. SU(3) breaking corrections, which are common in 
A(B+ -> J D°/^ + )/A(fi+ -> L>°7r+) and -> D°K S )/A{B° -> L>°7r°), cancel in Q 

and do not affect this estimate. 

It is more difficult to obtain reliable estimates for r c and r n . The two parameters are 
expected to be smaller than one since they contain a CKM factor | V u bV cs /V c bV us \ ~ 0.4. 
The parameter r c involves also an unknown color-suppression factor in b —>■ ucs, while r n 
involves the ratio of color-suppression factors in b — > ucs and b — > cus. Since the dynamics 
of B — > Z)i^ decays (caused by 6 — > cus) and £> — > .Dif (caused by 6 — ► ucs) are different, 
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color-suppression may be different in the two cases. This introduces large uncertainties in 
r c and r n j^ . 

It is easier to justify an approximate relation between the magnitudes of the two color- 
suppressed amplitudes A c r c and A n r n . Noting that the two processes B + — > D°K + and 
B° — > D°K° differ only by the flavor of the spectator quark, one expects 

A c r c ~ \/2A n r n . (5) 

In Section V we will discuss the approximation involved in this relation and a way of testing 
it experimentally. This approximate equality implies that the sensitivity to 7 is comparable 
in charged and neutral B decays, since the sensitivity in each case is governed by the smaller 
of the two interfering amplitudes. This point provides a major motivation for our study. 



III. TWO-BODY AND QUASI TWO-BODY D DECAYS 

Considering decays of D° and D° into a generic two-body or quasi two-body hadronic 
state fn and its CP conjugate fo, we denote the corresponding amplitudes by 

A(D° -> f D ) = A(D° -> f D ) = A f , 

A(D° -> f D ) = A(D° -> f D ) = A f r f e iS f , (6) 

where by convention Af > 0, r/ > 0, and < 5f < 2ir. Here and below we set the weak 
phase in D decays to zero and neglect D° — D° mixing. The effects of D° — D° mixing can 
be included as in j^j, but are not further discussed here. 

Using these notations, one finds expressions for decay rates in charged B decays, 

T(B + - f D K + ) = A\A) [1 + r 2 c rj + 2r c r f cos(S c + 5 f + 7)] , 

r(£T -> f D K-) = A\A) [1 + r\r) + 2r c r f cos(5 c + 5 f - 7)] , 

T(B + - f D K + ) = A\A) [r 2 + r ) + 2r c r f cos(5 c - 5 f + 7)] , 

T(B- -> /dA-) = A c 2 4 [r 2 + rj + 2r c r/ cos(5 c - 5/ - 7)] • (7) 

Combining B + and -B - decay rates for states involving a common D decay mode, fo or f^, 
one finds 

(r(B - / D ^c)) = r(5+ -> f D K+) + T(B- - fjjK~) 

= A\A) [(l + r c 2 )(l + rJ)+4r c r / cos(5 / + 7)cos5 c ] , 

(t(b - ; D ^ C )> = r(fi + - fr,K + ) + r(s- - ; D in 

= A 2 AJ [(l+r c 2 )(l + rJ)+4r c r / cos(5 / -7)cos5 c ] • (8) 
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Studying neutral B decays, one finds similar expressions for untagged decay rates 



0: 



(T(B - f D K n )) = T(B° -> JdKs) + r(5° -> fjjKs) 

= A 2 n A) [(l+r2)(l + rJ)+4r n r / cos(5 / + 7)cos(5 n ] , 

(T(B -> /^ n )) = r(5 -> JdKs) + r(5° -> fnKs) 

= A 2 n A) [(1+0(1 + rj) + 4r n r/C0s(5/- 7 ) cos 5 n ] . (9) 

Individual time-dependent decay rates for -B°(t) — > fnKs, B°(t) — > /dKs and their CP- 



conjugates are given in |20j, and include more information than the untagged rates. These, 
however, will not be needed in the following. 

The decay rates in Eqs. (jSJ) and (J2J) display a dependence on two types of quantities. Am- 
plitudes and strong phases in B — ► Dif, (Aj, r i; z = c, n), which in general obtain different 
values in charged and neutral B decays, and the corresponding quantities in D°/D° — > f D , 
(Af,rf,5f), which are common to both B + and B° decays. We will refer to these quantities 
as B and D decay parameters, respectively. In the following we will assume that the D decay 
quantities Af and Tf have been measured and are known. They can be obtained through 
branching ratio measurements in an independent sample of neutral D mesons, flavor-tagged 
through their production in the decay D* + — ► D n + [30]. In quasi two-body D decays the 
three parameters Af,Tf and Sf can be determined simultaneously thr oug h a complete Dalitz 
plot analysis. Although the phases 5f can in 

principle be measured 0H HQ 

, we will 

treat them as unknown, unless indicated otherwise. 

Note that the combined B ± rates in (jHJ) and those in the untagged B° decays (JHJ) depend 
in each case only on two combinations of B decay parameters, 

Xi = A 2 {I + r?) , Y t = 2A 2 n cos £ , i = c,n, (10) 

which obey 

\Yl\<Xi. (11) 

We see that the individual branching ratios for B° — > D°Ks and B° — > D°Ks, proportional 
to A 2 n and A^r 2 , respectively, cannot be measured from untagged decays alone. 

The three D decay parameters Af,rf and 5f in (JHJ depend, of course, on the final state 
fu. One may distinguish between three cases for which we give examples: 

1. f D = CP-state (e.g., f C p+ = K + K~, / CP _ = K s n°), for which r C p± = l,cos5 C p± = 
±1, 

2. /d = flavorless (e.g., K* + K~), for which r/ = 0(1) but generally r/ ^ 1,5/ = 
unknown, 
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3. fz> = flavor state (e.g., K + ir ), for which rf ~ tan 2 9c |3Q|, 5/ = unknown, where 9c 
is the Cabibbo angle. 

Using Eqs. and the observation (JTUJ), it is simple to show that 7 may be determined 
solely from untagged B° decays. Consider N different non-CP neutral D decay modes 
(k = 1, ...,N) together with their CP conjugates fj^, as the final states in the -B — > DK$ 
decay chain. The unknown variables are 7, X n , Y n and N strong phases 8h. Eqs. (JHJ), which 
provide 2N measurables for iV + 3 unknowns, are solvable for N > 3. That is, 7 may be 
determined from untagged B° — > DK$ decay rates, where D° is observed in at least three 
different non-CP decay modes and their CP conjugates. This argument may be generalized 
to include other untagged B° decays, such as B° — > D*K$ (D* — > D°n°). Assuming M 
different B° decay modes of this kind, each of which introduces a pair of unknowns X 3 n and 
Y£ (j = 1, M), one has 2MN measurables for 2M + N + 1 unknowns. For M > 2 this set 
of equations is solvable for N > 2. Namely, two non-CP decay modes of D° are sufficient 
for determining 7 from untagged B° — > DKs and B° D*K$. 

For a CP-eigenstate the strong phase 5/(cp) is either or n. In this case the two equations 
in (jHJ) become identical and provide a single measurable. Choosing the decay modes to 
be (i) an even-CP state, (ii) an odd-CP state, and (iii) a single non-CP eigenstate and its CP 
conjugate (involving an unknown phase 5f), one can solve the four equations for 7, X n , Y n 
and 5f. For this case we now derive an explicit expression for tan 2 7 in terms of measurable 
rates. The derivation holds for both charged and neutral B decays. 

Using Eqs. (jSJ) and (0), one has 

(T(B - fcp±Ki)) = 2A 2 CP± [X t ± Yi cos 7] , (i = c, n) , (12) 

where the two signs on the right-hand-side correspond to positive and negative CP- 
eigenstates. Adding and subtracting rates for even-CP and odd-CP eigenmodes, one finds 

, ^ (T(B — > fcp+Kj)} (T(B -> fcP-Kj)) 

^ s m^ ± m.m^=m = , yimi , (i= , B) . (14) 

Z/1 CP+ Z/1 CP- 

These definitions apply in practice to sums over individual CP states. Eq. (jT3j) provides 
the most direct way to determine X c and X n . We define CP- conserving rate asymmetries 
between even and odd CP-states, 

4, p = |p, (i = c ,n). (15) 

^CP 

Using Eq. one obtains two potential inequalities for cos 7 in terms of these ratios, 

|cos7| > I-4-cpI > (i = c,n) . (16) 
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This inequality holds se par ately for charged and neutral B decays. The inequality for neutral 
B decays was noted in [22]. 

As mentioned above, in order to determine 7 from untagged neutral B decays, one needs 
in addition to the two rates for CP eigenstates given in (|T2"|) (i = n) two rate measurements 
for a non-CP state fr> and its CP-conjugate Jd- These rates are given by the two equations 
in ©. We denote the sum and difference of these rates and of the corresponding rates in 
charged B decays by 

yi _ (T(B -> f D Kj)) + (T(B -> f D Kj)) 
f ~ ^/( 1 + r /) 

^ ""^1!.;^" ' li=c - nh (17) 

which imply CP-violating asymmetries, 

4 = A}/E}, (z = c,n). (18) 

It is then straightforward to show that tan 2 7 is proportional to {A l j) 2 and is given by 

t 2 (4) 2 _ (4) 2 (S>/Sbp) 2 v nq , 

tan7 p 2 (A^) 2 - (S^ P - p 2 (^cp) 2 -(1-^p) 2 ' l * C ' nJ ' UJj 

where 

" s rr% • (20) 

This result applies to both charged and neutral mesons. We stress that when determining 
7 we do not rely on separating the two terms A 2 and A 2 rf contributing to X{. 

While the possibility of measuring 7 from untagged B° decays alone is interesting, the 
most efficient way to determine the weak phase would be to combine information from decays 
of charged B decays and untagged neutral B decays. The derivation we have just presented 
applies also to charged B decays alone. One needs to measure only the combined B^ rates 
given in Eqs. (JSj). without needing to separate the small A\r 2 c term. Thus, Eq. (|T9|) also 
gives tan 2 7 in terms of decay rates for combined B^ — > DK^ events, where D mesons are 
observed in decays into an even-CP, an odd-CP eigenstates, and a non-CP flavorless state 
or a flavor state. We see that, in principle, a determination of 7 does not require measuring 
CP asymmetries in B ± — > fcpK ± . These asymmetry measurements, for even-CP and odd- 
CP states [3 5 provide additional useful information. Using all these measurements 
together with rate measurements for untagged neutral B decays will lead to a more accurate 
determination of 7 than when using only charged B mesons. 

Eq. (|19j) displays an explicit dependence of 7 on rate measurements defined in 
Eqs. (fT3j) . ([H)) and (JT7J). We see that tan 7 is proportional to the CP asymmetries 
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Af ~ Tf cos Si sin Sf sin 7 (z = c, n), indicating that the sensitivity for measuring 7 increases 
with rj. The sensitivity depends also on the value of Sf. In the extreme case that Sf vanishes 
(mod 7r) the two rates in (jSJ) (and in Q) become equal and 7 cannot be extracted. For 
the two-body flavor state / = K + tt~, the phase Sf vanishes in the SU(3) symmetry limit. 
However, SU(3) breaking effects are known to be large in D decays. Con sequ ently sizable 
values of Sf have been calculated for this final state in several models |34|. This phase 
can be measured at a charm factory 0,E3|- There are experimental indications for small 
phases in two cases of quasi two-body states, / = K* + n~ where Sf = (12 ± 3)°(mod it) was 
measured and / = p + n~ where Sf = (4 ±3 ±4)° was measured [35]. These results were 
obtained by studying the Dalitz plot of D° — > Ks^ + ^ and D° — > 7r + 7r~7r°, respectively. No 
measurement exists for Sf in D° — > K* + K~ , which can be measured by studying the Dalitz 
plot of D° -> K + K-7i° H- 

As we note in the next section, a complete Dalitz plot analysis of three body D° decays 
involves other strong phases, which are large in regions where two resonances overlap or 
when resonances interfere with non-resonant contributions. These phases will be shown to 
be useful when studying 7 in B — > DK, where the D meson is observed in a three-body 
final state. 



IV. B -> DK S OBSERVED IN MULTI-BODY D DECAYS 

The study of untagged B° — > DK$ presented in the previous section for D mesons 
decaying in two-body modes may be extended to multibody decays. To be specific, we focus 
on the case of the three-body D decay, 

D -> K s n-n + , (21) 

n 

following the discussion of B ± — > DK ± in [9J. In order to make our point, we start with 
a model-independent approach. We also explain how modeling the amplitude for D° — > 
Ks7r~7r + in terms of a sum of a given set of intermediate resonances [3^, as done recently 
in may help in reducing the experimental error in 7. 

We denote the amplitude for D° — > K s ti~'k + at a given point in the Dalitz plot by 

A(D° - K s ( Pl )n-(p 2 )n + (p 3 )) = A(s 12 , Sl3 )e i<5(si2 ' si3) , (22) 

where = (pi +pj) 2 . As in two body decays, we use the convention A(si2,sis) > and 
< S(s\2, S13) — 27T. Assuming that CP is conserved in this decay, one has 

A(D° ^ K s ( Pl )7r-(p 2 )n + (p 3 )) = A(D° ^ K s ( Pl )n-(p 3 )n + (p 2 )) = A(s 13 , s 12 ) e l5 ^ s ^ . 

(23) 



9 



That is, the (complex) decay amplitude for D° at a given point (si2,si3) in the Dalitz 
plot equals the decay amplitude for D° at a point (S13, Si 2 ) obtained by reflection across a 
symmetry axis corresponding to exchanging the momenta of the two pions. 

The density of events in the D decay Dalitz plot for untagged B — > (Ks^^ + )dKs is 
obtained using Eqs. (0), (|2*2*j) and (or from the time- dependence in Appendix A), 

d 2 r 



[B°/B° - [/^bi)7r-(p 2 )7r + (p 3 )]D^) = ^ [(A 2 (s 12 , sis) + A 2 ( Sl3 , (1 + 



ds\ 2 ds\ 3 

+ 4r n A(si 2 , si 3 )v4(si 3 , Si 2 ) cos(5(si 2 , s 23 ) - £(si 3 , s 12 ) + 7) cosS n ] . (24) 

This is in complete analogy with the first of Eqs. ©• The second equation, describing 
the density at the point of reflection across the symmetry axis, involves an opposite sign 
for 7. Integrating ()24j) over an area (a bin) i lying below the symmetry axis and over a 
corresponding symmetry-reflected area % lying above the symmetry axis, one has 1 

(r f > = [dT(B°/B - [K S 7i-7i + ] D K s ) = X n (Ti + Tj) + 2Y n [ Cl cos 7 - s< sin 7] , 

J i 

(T-)= [dT(B°/B° ^[K s n-n + } D K s )=X n (T i + T\) + 2Y n [c l cos 1 + s i sm 7 } , (25) 

J i 

where we define 

Ti = ds 12 ds 13 A 2 (si 3j s 2 3) , 



ds 12 ds 13 A(s 12 , s 13 )A(s 13 , s 12 ) cos(5(si 2 , s 23 ) - 5(s 13 , s 12 )) , 

Si = J^ds 12 ds 13 A(s 12 , s 13 )A(s 13 , s 12 ) sin(5(si 2 , s 23 ) - S(s 13 , s 12 )) . (26) 

The partial rates Tj in D decays may be measured using flavor-tagged D° decays and are 
assumed to be known. The other D decay variables, q and Sj, which in principle can 
be measured model-independently at a charm factory (up to a sign ambiguity in Sj), will 
nonetheless be taken as unknown. Consider k different bins % lying below the symmetry 
axis, each contributing two unknowns and Sj. Together with X n , Y n and 7, there are 
2k + 3 unknowns. Eqs. (|25|). which provide 2fc measurables ((Ti) and are therefore 

unsolvable. 

The situation changes when one measures another neutral B decay of this type, e.g. 
the sequence B° — > D*K$, D* — > D°ti , D° — ► Ks^ + ^ ', which introduces a pair of new 
variables analogous to A ra , y n . In this case one has 4/c measurables for 2/c + 5 unknowns, a 



1 The method outlined here applies to any multibody D° decay with the set of equations l|25(l unchanged. 
If the i-th bin is in the phase space of a final state fu, then the z-th bin is a CP transformed bin in the 
phase space of the corresponding fjj state. 
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solution for which requires k > 3. That is, 7 may be determined by measuring partial rates 
in the two untagged neutral B decay modes, B — > DK$ and B — > D*Ks, for at least three 
pairs of Dalitz plot bins in D — > fTs7r + 7r~ . 

A more powerful approach is to combine information from all the untagged neutral i? de- 
cays and charged B decays, using both multibody and two-body D decays. For instance, in 
analogy with Eqs. (J7J), the decays B ± — > {Ksti + ^)dK ± provide four measurables for each 
bin ji.9], instead of the two in (|2*5|l. Combining charged and untagged neutral B — > DK de- 
cays, where D — > fCs7r + 7r~ , yields 6/c measurables for 2/c + 6 unknowns, q, Sj, v4 c , r c , <5 C , X n , F„ 
and 7. Therefore, two pairs of bins provide an overconstrained system of equations for de- 
termining 7. 

The two equations (|25|) are not mutually independent when = 0, in analogy with 
the singularity noted in Eq. (|19|) when 5f = (mod tt). However, the relevant strong 
phase differences which determine Sj are large at least in some regions of the Dalitz plot 
of D° — > Ks'k+'k^ . Consider, for instance, the two overlapping regions of a vertical band 
describing the Cabbibbo-allowed mode K*~tc + and a horizontal band describing the doubly 
Cabibbo-suppressed mode K* + tt~ with a diagonal band representing the Cabibbo-allowed 
mode p°K°. The local strong phases which determine for these two regions are the phase 
differences between amplitudes describing the sum of the K*~ir + and p°K° contributions 
and the sum of the K* + tc~ and p°K° contributions. These phases are large and vary a lot 
over the overlapping regions because of the two largely different K*it contributions in the 
two amplitudes. 

One can reduce the number of unknowns appearing in the determination of 7, if the un- 
knowns coming from the D decay, Cj and s,, appearing in ()25|) . are determined independently. 
This can be done by assuming a Breit-Wigner (BW) form for the intermediate resonances 
contributing to this decay 9]. The parameters of the model describing the D decay ampli- 
tude can then be fitted to data of tagged D decays, which are abundant at 5-factories. The 
observables in ()25|) now depend only on three unknowns, X n , Y n and 7. 

It is hard to quantify the theoretical error introduced by assuming a BW form. One way 
to proceed is to change the number of resonances and see how the sensitivity changes. This is 
only a partial determination of the error. Another source of error is the accuracy of the BW 
assumption. This can be determined by the goodness of the fit to the tagged D decays, or by 
using a different mode, for resonances, such as a Kuuatrix ,uode. for wide resonances ft. A 
rough estimate of the theoretical error caused by assuming a superposition of BW amplitudes 
is about 10° Further studies are required in order to evaluate possible contributions of 
non-BW terms in the D decay amplitude and their effect on determining 7. As mentioned, 
this model-dependence can be avoided by measuring the parameters q and Sj at a charm 
factory. 
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V. USING ISOSPIN AND NEGLECTING AN ANNIHILATION AMPLITUDE 



As already mentioned, the most powerful approach for improving the determination of 
7 is to combine charged B decays with the information from untagged neutral B decays. 
Adding untagged neutral B decays to a sample of charged B decays with the same D 
decay final states introduces only two unknown parameters, X n and Y n . Here we discuss 
an approximation which may be used to reduce the number of parameters further. This 
introduces a theoretical error in 7. Yet, it is worthwhile considering such an approximation 
as long as this error is smaller than the statistical error. 

We recall two isosprn relations 0, one for 5 a. transitions, 

A(B° -> D~K + ) = A(B + -> D°K+) - A(B° -> D°K°) , (27) 

and another for b — > ucs transitions, 

A(B° -> D°K°) = A(B + -> D°K + ) + A(B + -> D + K°) . (28) 

The amplitude of B + — > D + K° is pure annihilation, and is expected to be smaller than the 
other two amplitudes in the last relation j^J. The absence of rescattering effects, which may 
enhance this amplitude to a level comparable to the other two amplitudes in this relation, 
can be tested [39j] by setting very stringent experimental bounds on the branching ratio for 
B + -> D + K°. Neglecting A(B+ -> D + K°) } Eq. fl2EJ) reduces to Q 

A(B° -> D°K°) = A(B + -> D°K + ) . (29) 

Namely, the two color-suppressed amplitudes have equal magnitudes and equal strong 
phases. Note that the error due to isospin violation, at most a few percent, is likely to 
be much smaller than that involved in neglecting the annihilation amplitude. 

Equations ()27|) and (}2*9"j) may be used to simplify the determination of 7 when combining 
B^ 1 decays and untagged B° decays. The two equations imply 

A\ + 2Al - 2V2A c A n cos(5 n - 5 C ) = T(B° -> D~K+) , 

V2A n r n = A c r c . (30) 

The right-hand side of the first equation, which involves a color-allowed process, has already 
been measured 0| and will be assumed to be given. Eqs. (}30|) reduce the six parameters 
describing charged and neutral B — > DK decays, A^r^d^i = c,n), to four independent 
ones. The measurable parameters in untagged B° — > DK S decays, X n = A^(l + r^) and 
Y n = 2A 2 n r n cos5 n , can now be expressed in terms of the three B + decay parameters and a 
single B° decay parameter (8 n , for instance). That is, under the above assumption, adding 
information from untagged neutral B° decays to studies of B ± decays involves a single new 
unknown parameter instead of two parameters. This is expected to reduce the statistical 
error in determining 7. 
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VI. CONCLUDING REMARKS 



Before concluding let us make several general comments: 

• Since the data set of B + — > DK + is larger than that of B° — > DK$, hadronic param- 
eters such as X c are easier to measure than X n . As noted, the approximate relation 
\p2A n r n ~ A c r c implies comparable sensitivities to 7 in charged and neutral B decays. 
The sensitivity in B° decays is smaller by a factor v^2 since only K$ mesons are ex- 
perimentally useful. Moreover, the detection efficiency for K$ is about a factor of two 
smaller than that for charged kaons. This is expected to reduce somewhat the effect 
of neutral B decays on determining 7. 

• Our study focused on B —>■ DK decays. It can be extended to multibody B decays, 
including B + — > D*K + and B° — > D*K S , where D* — > Dn°, as well as to the self- 
tagged decays B + -> D^K*+ and B° -> D^K*°. This would add to the statistical 
power of the analysis since the parameters related to D decays are common to these 
processes and to B — > DK. 

• The method we discussed for B° decays can be applied also to B s decays, replacing 
the K s by <fi, rj', 77. In that case the advantage of being able to use untagged data 
is greater, because the hadronic environment where B s decays will be studied makes 
flavor tagging less efficient. We have neglected the width difference between the two 
neutral B meson states, which is a very good approximation for nonstrange B mesons. 
In the case of B s , the width difference is expected to be nonnegligible and may be 
taken into account in a straightforward manner. 

• In our discussion we assumed that strong phases in D decays are unknown and need 
to be determined from the analysis simultaneously with 7. As we did already men- 
tion, the stron g ph ases in two-body and quasi two-body D decays can be determined 
independently \WL 31, 32, The strong phases in three-body decays may also be 



determined by assuming that the D decay amplitude is given as a sum of Breit-Wigner 
amplitudes. Knowledge of strong phases would imply, for instance, that fewer D decay 
modes are needed in order to determine 7 from untagged decays alone. In practice, 
this implies that a combined fit of the data to fewer hadronic parameters will result 
in a smaller error in 7. 

The extraction of 7 from B ± — > DK^ involves a number of discrete ambiguities 0, 
l^lj l . The ambiguities in untagged B° decays may be identified in Eqs. (jOJ) which are 
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invariant under 



P ox = {7 -> 5 f , 5/^7}, 

P_ EE {7 -> -7 , 5f -> -<5/} , 

-Ptt = {7 7 + TT , 5/ -> 5/ + 7T or 5 n -> S n + 7r} . (31) 

Once several two-body D decay modes are combined, or once a multibody D decay 
mode is used, the first two ambiguities may be resolved. The ambiguity P ex is lifted, 
since 5/ is not expected to be the same for all two-body and quasi two-body D decay 
modes, and is known to change over the Dalitz plot in three-body decays. Measuring 
the sign of Sj in (|25|) through a fit to a sum of Breit-Wigner resonance functions 
would resolve the P_ ambiguity. This introduced essentially no model-dependence, 
since one needs only the sign of Sj, which is easily determined in the vicinity of a BW 
resonance. Resolving P_ avoids the ambiguity 7 — > 71 — 7, a combination of P_ and 



P^, which is particularly problematic in view of the proximity of 7 to ir/2 [4l|. The 
only remaining ambiguity is 7 — > 7 + 7r. This ambiguity is the least problematic, since 
the two corresponding values of 7 are maximally separated. 

To conclude, we have studied the information obtained from untagged neutral B decays 
of the type B — > DK$ involving several D decay modes. We have shown that these mea- 
surements alone can, in principle, determine 7. Of course, B° tagging information, while 
limited, can only improve this determination. By combining information from untagged B° 
decays with that obtained in corresponding charged B decays one gains statistics, thereby 
permitting a more accurate determination of 7. While statistics are limited, one may neglect 
an annihilation amplitude in B — > DK , reducing by one the number of hadronic parameters 
and resulting in a smaller experimental error in 7. This introduces a theoretical error in 7 
which must be further studied. 
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APPENDIX A: TIME-DEPENDENT B -» DK S WITH MULTIBODY D DECAYS 

In this appendix we provide a formalism allowing the extraction of 7 from time- dependent 
rates in B° —>■ JdKs where fo is a multibody final state. As we show, this does not only 
serve the purpose of determining 7, but also helps resolve the current two-fold ambiguity, 
(3 — > 7r/2 — (3. For simplicity we take fo to be the three-body final state J\~s7r + 7r~ studied 
in Section IV. 

Time-dependent partial rates, integrated over a bin i in the Dalitz plot of D — > J<s7r + 7r~ 
lying below the symmetry axis, and over a corresponding symmetry-reflected bin % above the 
axis, are readily calculated for initial B° and B° states (a positive B meson bag parameter 
is assumed 

= [dT (B°(t) -> (K s n-n + ) D K s ) = e~ FBt A 2 n x 

J i 

{/+cos 2 + Jrsin 2 f^^j +S lS m(Am B t)} , (Al) 

T t = [dT (B°(t) - (K s 7f-n + ) D K s ) = e~ VBt A 2 n x 

J i 

[i-; cos 2 (^y^) + 7+ sin 2 (^^) - 51 sin( Am B i) } , (A2) 
T-i = [dT (B°(t) -> (^ s 7T-7r + ) D /i" s ) = e- rsi A 2 x 

I cos 2 + ir sin 2 + Sj sin(Am fl t) } , (A3) 

r 5 = [dT (B°(t) -> (^7r-7r + ) D ir s ) = e" rs M 2 x 

{/- cos 2 + I-t sin 2 - S 5 sin( Am B t) } . (A4) 

The six observables determined from the time-dependence, I i , L , Si and Sj, are defined in 
terms of the quantities in Eqs. 



if = Ti + r 2 Ti + 2r n [cos(7 ± 5 n )ci =F sin(7 ± 5 n )si] , (A5) 
$ = r n T- sin(2/5 + 7 - 5 n ) + [sin(2/5) Ci - cos(2/5) Si ] 

+ r 2 [sin(2/3 + 2 7 )q + cos(2/3 + 2 7 ) Si ] + r n T, sin(2/3 + 7 + 5 n ) . (A6) 

Expressions for the observables L and Sj are obtained from (jA5|) and ()A6|) by replacing 
Ti <-> Tj, Sj — > — Sj. Note that J,j correspond directly to the partial decay widths r.-= for 
B ± -> fJ^ defined in 0. 

Dividing the Dalitz plot into k pairs of bins, i and i, the 6k observables permit an 
extraction of 7. There are 2k + 4 unknowns, Cj, Sj, A n , r n , S n , 7 (sin(2/5) is assumed to be 
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known), so that the system is solvable for k > 1. Namely, in order to determine 7 from 
time-dependent decay rates into {Ks'k + t^~)dKsi it is sufficient to divide the D decay Dalitz 
plot into two bins, symmetric with respect to the symmetry axis. 

The solution for 7 involves a four-fold discrete ambiguity. Equations (jAl|) - (jA4|) are 
invariant under the following four independent discrete transformations: 

PJ = {7 -»• 7 + S n -> 5 n + tt} , 

= {7 -> 7 + 7T, /3 — ► /3 + 7T/2, Cj -> -C;, Sj -> -Si} , 

P_ = {7 -> -7, /? tt/2 - (3, 5 n -> -<5 n , s< -> -s,} . (A7) 

The P^ ambiguity can be resolved model-independently, either by using the sign of sin 2(3 
or by measuring the sign of q at a ^(3770) charm factory [9|. The P_ ambiguity can 
be resolved if one determines the sign of Sj by fitting the Dalitz plot to a sum of Breit- 
Wigner forms. (See Section VI.) Note that resolving the P_ ambiguity in this way leads to 
the determination of the sign of cos (2f3) in an essentially model-independent way. It also 
determines the sign of 7 or equivalently the sign of cos(2a). Fixing the sign of cos(2/3) is a 
consequence of knowing the sign of in multibody decays. This is impossible in two-body 
D decays, where the sign of sin 8f cannot be determined. The remaining two ambiguities, 
F*2 and Pf , cannot be resolved without further theoretical input. 
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